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Abstract 

We consider Diophantine quintuples (a, b, c, d, e}, sets of integers with a<b<c<d< 
e the product of any two elements of which is one less than a perfect square. Triples of 
the first kind are sets {A, B , C} with C > B 5 . We show that there are no Diophantine 
quintuples {a, 6, c, d, e } such that (a, b , d} is a triple of the first kind. 


1 Introduction 

Define a Diophantine m-tuple as a set of m distinct positive integers cq < cq < ■ ■ ■ < a m 

such that GqOj + l is a perfect square for all 1 < i < j < m. For example, the set {1, 3, 8,120} 

is a Diophantine quadruple. Throughout the rest of this article we simply refer to m-tuples, 
and not to Diophantine m-tuples. 

One may extend any triple {a, 6, c} to a quadruple {a, 6, c, d + } where 

d + = a + b + c + 2 abc + 2 rst, r = V ab + 1, s = y/ac + 1, t = Vbc + 1, 
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by appealing to a result by Arkin, Hoggatt and Straus [2]. Such a quadruple {a,b,c,d + } is 
called a regular quadruple. Arkin, Hoggatt and Straus conjectured that all quadruples are 
regular. Note that any possible quintuple {a, b, c, cl, e} contains the quadruples {a, b, c, d} 
and {a, b, c, e}. If the conjecture by Arkin, Hoggatt and Straus were true then d + = d = e, 
whence d and e are not distinct. This implies that there are no quintuples. A partial result 
towards this is the following theorem, proved by Fujita HD], 

Theorem 1 (Fujita). If {a,b,c,d,e} is a Diophantine quintuple with a<b<c<d<e, 
then d = d + . 

When attempting to bound the number of possible quintuples {a, b , c, d, e} it is useful 
to examine the relative size of b and d. To this end, Fujita mi considered three classes of 
triples {A, B,C}, namely, triples of the first kind in which C > B 5 ; triples of the second 
kind in which B > 4 A and B 2 < C < B 5 ] and triples of the third kind in which B > 12 A 
and B 5 / 3 <C <B 2 . 

Lemma 4.2 in [6j states that every quadruple contains a triple of the first, second, or 
third kind. Specifically, we have 

Lemma 1 (Lemma 4.2 in [6]). Let {a,b,c,d,e} be a Diophantine quintuple with a < b < 
c < d < e. Then 

1. {a, b, d} is a triple of the first kind, or 

2. (i) {a, b, d} is of the second kind, with Aab + a + b < c < b 3 ^ 2 , or 

(ii) {a, b, d } is of the second kind, with c = a + b + 2r, or 

(in) {a, b, d} is of the second kind, with c > b 3 ^ 2 , or 

(iv) {a, c, d} is of the second kind, with b < 4a and c = a + b + 2 r, or 

3. {a, c, d} is of the third kind, with b < 4a and c = (Aab + 2) (a + b — 2 r) + 2(a + b). 

Lemma [Tj allows us to take aim at various triples: to prove there are no quintuples one 
need only prove the nonexistence of quintuples containing each of the kinds of triples listed 
in Lemma [0 The aim of this paper is to show that there are no quintuples satisfying part 1 
of Lemma [0 

The current best estimate^], by Elsholtz, Filipin and Fujita [6], is that there are at most 
6.8 • 10 32 quintuples comprising 

• 5.05 • 10 15 possible quintuples derived from triples of the first kind, 

• 6.72 • 10 32 possible quintuples derived from triples of the second kind, and 

• 1.92 • 10 26 possible quintuples derived from triples of the third kind. 

lr The second author has announced m that there are at most 2.1 • 10 29 quintuples, and that there are 

no quintuples of the kind 3 in Lemma 2] 
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Attempts to bound the total number of quintuples have used a result by Matveev [23 
on linear forms of logarithms. This is not required for triples of the first kind, which greatly 
simplifies the exposition. Rather, the bounds on b and d come from work of Filipin and Fujita 
[7] wherein so-called gap principles between solutions of Pell’s equation are considered. 

In £j2]we improve Fujita and Filipin’s proof slightly. This enables us to perform compu¬ 
tations in §[3]that prove: 

Theorem 2. There are no Diophantine quintuples {a, b, c, d, e} such that {a,b,d} is a Dio- 
phantine triple with d > b 5 . 

Independently of this work, Cipu [5] has considered the same problem. Indeed, his 
Theorem 1.1 states that “the quadruple left after removing the largest entry of a Diophantine 
quintuple contains no standard triple of the first kind”. This is somewhat stronger than our 
Theorem [2l though either result may be used to eliminate the triples in part 1 of Lemma [0 
Moreover, Cipu shows [Sj Thm. 1.1] that there are at most 10 31 quintuples. 

We were only made aware of this result when this paper was essentially completed. 
Elements of Cipu’s paper, such as the reliance on inequalities of the form (J2]) and (J7]) are 
similar to ours. Our computational approach is different, and, in particular, our Algorithm 
[2] introduces some new ideas to the field. 
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2 Bounds considered by Filipin and Fujita 

If a double or a triple cannot be extended to form a quintuple, we need not consider them 
in what follows. We call such doubles and triples discards. Many discards are known; we 
require only a few. The first, due to Fujita [9] (see also 0), that shows that the double 
{k, k + 2} is a discard. In addition, Filipin, Fujita and Togbe [8J Cor. 1.6, 1.9] proved that 
the following are discards for k > 1 

{3k 2 — 2k, 3£r + + 1}, {2(/c + l) 2 — 2(k + 1), 2{k + l) 2 + 2{k + 1)}, 

{{k + l) 2 — 1, (fc + l) 2 + 2{k + 1)}, {k,4k + 4}. U 

For the following arguments we make frequent reference to Lemma 2.4 and (2.9) in [7]. 
These rely on a property denoted as ‘Assumption 2.2’ concerning the relations between the 
indexed solutions of the Pcllian equations associated with the hypothetical quintuple. For 
ease of exposition, we do not write out the details of Assumption 2.2; we merely note that, 
for its application to quintuples containing triples of the kind in 1 in Lemma [H it is satisfied 
as per [3, p. 303]. 

We also note, for the convenience of the reader, that a' in [7, (2.9)] is defined as a' = 
max{6 — a, a}. We assume that {a, b, d} is a triple with d > b 5 . We consider two cases 
according as b > 2a or b < 2a. 
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2.1 When b >2a 

Combining (2.9) and Lemma 2.4(i) in [7] we have 


0.178a 1//2 d 1//2 log(4.001a 1/,2 (6 — a) 1/,2 6 2 d) log(1.299a 1//2 6 1/,2 (6 — a) 1 d ) 

46 < log(46<f) log(0.1053a6 _1 (6 — a)~ 3 d) 


( 2 ) 


It is easy to see that both factors of the numerator and the second factor of the denominator 
are increasing in a for b > 2a. Hence, given that 1 < a < 6/2, the right side of (J2J) is bounded 
above by 

log(2.00056 3 d) log(1.8371d) 
log(46d) log(0.10536 _4 d) 


It is easy to see that (logH 1 x)(logH 2 x)/((logH 3 x)(logH 4 a;)) is decreasing in x if A 4 > A 3 
and A 2 > A 4 . Since 6 > 2 it follows that (J3J) is decreasing in d, and since d > 6 5 , we have 
from ([2]) that 

b s /2 4 log(2.00056 8 ) log(1.83 716 5 ) 

0.178 log(46 6 )log(0.10536) ' U 

We find that flD) holds provided that 6 < 50; in [7] the bound derived is 6 < 52. Thus we 
need only consider those pairs {a, 6} for which 1 < a < 6/2 < 25. We enumerate these pairs, 
and test them against the inequality in ([2]). After discarding pairs containing {k, k + 2} and 
those doubles in (0Q) we find that the only possibilities are 


{1,15}, {1,24}, {1,35}, {2,24}, {3,21}. (5) 

For each of these five doubles we insert values of a, b into (l2]l . and solve for d. For example, 
with {1,15} we find that d < 5.2 • 10 6 . We now search for all those d with 15 s < d < 5.2 • 10 6 
such that {l,15,d} is a triple. We find only one such value of d, namely d = 2030624. 
We now search for all those c G (15,2030624) such that {1,15, c, 2030624} is a quadruple. 
This yields exactly one value of c, namely c = 32760. Thus, one possible quadruple is 
{1,15,37260,2030624}. We continue in this way with each of the doubles in ([5]). We find 
that the only possible quadruples are 


{1,15,37260,2030624}, {1,24,148995,14600040}, {1,35,494208,70174128}. (6) 

Note that each of the quadruples in (0 has a — 1. The proof of Lemma 2.4(i) in [7], 
that leads to the ‘0.178’ on the right side of ([2]) can be improved significantly if it is known 
that a — 1. Indeed, if one there assumes n < 0.456 _1 c 1 / 2 one obtains (2.4) in [7] as well 
as the desired contradiction, as before. Thus, for a = 1 one may replace the ‘0.178’ in ([2]) 
by ‘0.45’. It now follows that the double {1,15} when extended to a triple {1,15, d} must 
have d <2 ■ 10 6 , whence we eliminate the first quadruple in (JH]). The other two quadruples 
are similarly eliminated. We conclude that there are no quintuples with b > 2a for which 
{a, 6, d} is a triple of the first kind. 
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2.2 When b <2a 


In this case we combine (2.9) and Lemma 2.4(iii) in [Tj and obtain 

a -i/ 2 ^i /8 log(4.001afe 2 d) log(1.299a 1 / 2 6 1 / 2 (6 — a)~ l d ) 

4 < log(4M) log(0.1053d6 _1 (6 — a) -2 ) 


( 7 ) 


We need a lower bound on b — a: given Fujita’s result that {k, k + 2} cannot be extended to 
a quintuple, we can write b — a > 3. Since, again, the logarithms dependent upon a in (J7|) 
are increasing with a, and since a + 3 < b < 2a we rewrite (J7]) as 


6 i/s 1 log(4.001& 8 )log(0.4336 6 ) 
log(45 6 ) log(0.42125 2 ) 


We find that (JS]) is true provided that b < 4.69 • 10 9 . In fact, we can squeeze a little more 
out of the argument. Let a G [1/2,1) be a parameter to be chosen later. Then, for a < ab 
we have 


i/8 i /2 log(4.001afr 8 ) log(1.299a 1 / 2 6 6 (l - a) *) 

log(46 6 ) log(0.42126 2 ) 


(9) 


whereas, for a > ab we have 


6 1 / 8 < 4 


log(4.0016 8 ) log(0.4336 5 ) 
log(46 6 ) log(0.10536 2 (l — a) -2 ) 


( 10 ) 


Therefore 5 1 / 8 is less than the maximum of the right sides in (J9]) and (ITOji . We find that 
choosing a = 0.9862 gives b < 1.3 • 10 9 . Filipin and Fujita |7] proved that b < 10 10 : while 
our improvement is only slight, it makes the problem computationally tractable. 

Filipin and Fujita proved also that d < 5 9 ; this was improved in [6] to d < b 7 ' 7 . We use 
the weaker bound d < b 8 for computational convenience. 

We hope now to search for possible quadruples {a, b, c, d} with (a, 6, d} a triple of the 
first kind. We first enumerate all double {a, 6} with a > 1, a + 2<6<2a and b < 1.3 • 10 9 . 
For each such doubles we enumerate all c where c < b 8 such that {a, b, c} is a triple. We now 
appeal to Theorem [T] and compute d = d + . If b 5 < d < b 8 and a, b and d satisfy inequality 
©, then we add the quadruples to our initial list. 


2.3 Specific bounds for d 


There are a little under 11 million quadruples in our initial list (details are given in j|3]). We 
now propose a criterion against which to test these specimens. 

Each quadruple {a, 6, c, d} gives rise to a system of Pcllian equations the solutions to 
which are indexed by integers m and n — see [Tj pp. 294-295]. One obtains tighter bounds 
by showing that m and n must be of roughly the same size. Implicit in the proof of Lemma 
2.3 in [7] is the following problem. Given {a, b, d} a triple of the first kind, and given m,n 
with m > 3 and n > 2 we want to find good bounds on 


log(2.001 2 M) + 1; log(1.994a 1//2 6 1/<2 ) 
log(1.994 2 ad) 


w(a, 6, d , n), 


m 

— < 

n 
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( 11 ) 








say. Note that w is decreasing with n. Let n 0 be the smallest value of n such that (n 0 + 
1 ) Ah) < w(a, b, d, no) = 71 , say. 

Now let 72 be any number satisfying 72 > 7 ^, and let 73 be any number less than 


1 

7i 




1 72 ad + 1 

V ad -I - 1 



u(a,d,7i,7 2 ). 


( 12 ) 


Filipin and Fujita have 71 = 1.2 ,72 = 1 . 45,73 = 0.0033. This gives them a criterion against 
which to test quadruples provided that b > 1.45a. We should like to take 72 to be less than 
1.45, so that we can test more quadruples. We do this in the following lemma. 

Lemma 2. Assume that 72 a < b < 2a, and that {a, b , d} is a triple with d > b 5 . Then 


7 3 a 1/,2 & 1 d 1 / 2 log(4.001a6 2 d) log(1.299a 1 / 2 6 1//2 (6 — a) 1 d) . . 

4 < log(4M) log(0.1053d6 _1 (6 — a) -2 ) 

Proof. Using (1T2T) . and that fact that 72 > 7 i, the proof follows exactly the same lines as in 
the proof of Lemma 2.4(ii) in [Jj. We complete the proof by combining the bound on n with 
(2.9) in [7j. □ 


3 Computations 


We first present a simple result on triples. 

Lemma 3. Let a, b, r be positive integers with a < b such that ab + 1 = r 2 . That is, {a, b} 
is a double. Then all admissible c n such that {a, b , c n } is a triple are of the form 

„ % 2 n ~ 1 _ Vi ~ 1 

Cn 7 

a b 

where c n is an integer and the x n ,y n are integer solutions to 


bx 2 — ay 2 = b — a. 

Proof. For c n to be admissible, we have 


ac n = x 2 — 1 


and 

bc n = y 2 - 1 

for some integers x, y. We now simply eliminate c n . 


(14) 


□ 
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We now require an efficient algorithm to identify low-lying solutions to (fT4lh We start 
by dividing throughout by g = (a, b ) to get 

tfx 2 — cJy 2 = — a) 

and then write X = tfx, D = affif and N = b\b^ — cd) to get 

X 2 - Dy 2 = N. (15) 

Since D is not a square by construction, this is an example of Pell’s equation which has 
been widely studied (see, for example, [13]). To find solutions, we first use Lagrange’s PQa 
algorithm to find (u,v), the fundamental solution to 

X 2 - Dy 2 = 1. 

We then use brute forc^l to locate all the fundamental solutions to (fT5|) . Each such fun¬ 
damental solution potentially gives rise to two infinite sequences of solutions, one generated 
by 

X n +i + Vn+i^D = (. X n + y n \Z~b)(u + v'/D) 

and the other by 

X n +1 + Vn+i^D = (X n + y n VD){u - vVd). 

Thus we have a series of recurrence relations, indexed by i that will generate all possible 
solutions, which we denote 


•U,n+1 f Eiipijni Ui,n) and £/i,n+l f 2/i(*^i,n5 2 /i,n) 


for % E [0, / — 1]. 

We can now proceed as described in Algorithm [1] to find all quadruples {a, b, c, d} with 
0<a<b<c<d, a + 2 < b < 2a, b < 1.3 • 10 9 and b 5 < d < b 8 that satisfy inequality (J7J) . 

We note that iterating over the divisors of r 2 — 1 is more efficient than iterating over a 
and b and testing whether ab+1 is a perfect square. When factoring r 2 — 1 for this purpose, 
we first factor r ± 1 and merge the results. 

We implemented Algorithm [1] in Pari [3] (for the factoring) and in ‘C’ using GMP [[12] 
(for everything else). We ran it on 5 nodes of the Lhiiversity of Bristol Bluecrystal Phase 
III cluster [1] each of which comprises two 8 core Intel Xeon E5-2670 CPUs running at 2.60 
GHz. The total time used across these nodes was about 40 hours. 

There are 4, 038,480, 906 pairs {a, b} with ab + 1 a perfect square, a + 2 < b < 2a and 
b < 1.3 ■ 10 9 . From these pairs, we obtained 12,115,454,363 potential quadruples with 
b < c < d and b 5 < d < b 8 . Applying inequality flT]) eliminated all but 10,811,817 of these 
quadruples and these survivors formed our initial list. 

2 Theorem 6.2.5 of [T5] gives upper bounds on such a brute force approach. 


7 



for r <— 2 to 1.3 • 10 9 do 
w <— r 2 — 1; 
for a\w, a < r — 1 do 

b w/a; 

if b > 2a or b > 1.3 • 10 9 then continue Solve Pell’s equation 
bx 2 — ay 2 = b — a to get / base solutions; 
for i G [0, / — 1] do 

x «- Xj j0 ; 

V t— 2/i,oi 

t 2 -1 

c t- 2 — 

a ‘ 

while c < b do 

1 1- fxi(x,y)- 
V <- fVi{x,y ); 
x <— t; 

c t- 2 — 2 ; 

end 

d •<— a + fe + c + 2afcc + 2 ryJ ac + l-\/6c + 1; 
while d < b 5 do 

1 1 - fxi(x,y ); 

2 / fVi{x,y ); 

x <— t] 

— 1 

c «- 2 — 2 ; 

a ' _ 

d a + 5 + c + 2a&c + 2 r\J ac + l\/6c + 1; 

end 

while d < b 8 do 

if c e Z and a, 6, d satisfy inequality (|7J) then 
) output {a, 6, c, d} 

end 

1«- fxi(x,y ); 
y <- fyi(x,y)-, 
a: ^— t; 

c t- 2L — 1 


d ^ — ci -I - b -p c -p 2 abc 2vyJac ~|-1 \/ be T 


end 


end 


end 


end 


Algorithm 1: Producing the initial list. 



















We now apply Algorithm [2j based on inequalities (ed, m and m, to prune our 
initial list. To illustrate, consider the quadruple {a, b, c, d} = {8,15, 21736,10476753} which 
survives Algorithm [H We have 

31 

— > w(a , 5, d, 30) 

30 v , , , ; 


but 


32 

31 


< w(a, b, d, 31) 


1.0330... = 7i. 


Now 6 /a = 15/8 > 7 ^ so we will take 72 = 15/8. We now set 


73 = v(a, d, 71 , 72) = 0.32555 ... 


and we find that the left hand side of (El) is 49.67... and the right hand side is 2.852... so 
the inequality fails and we have eliminated this quadruple. 

When coded in Pari, Algorithm Stakes less than 20 minutes on a single core to determine 
that none of the quadruples in our initial list is admissible. This proves Theorem [2] 


for (a, b, c, d} G initial list do 

n 0 4— smallest n such that ^44 < w (a, 6 , d, rt); 

71 4 — w(a,b,d,n 0 y, 

72 4 — max(J,7?); 

73 4- v(a,d, 71,72); 

if a, b, d, 73 satisfy inequality (113} then output {a, b, c, d} 

end 

Algorithm 2: Pruning the initial list. 
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